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ABSTRACT 

We have obtained simultaneous optical and near infrared observations of 4U 
1543-47 during its 2002 outburst. The most striking feature of the outburst light 
curve is the secondary maximum which appears after the object transitions into 
the low-hard state. This secondary maximum is much stronger in the infrared 
bands than optical. We suggest that the origin of the secondary maximum flux 
may be synchrotron radiation associated with a jet. Close infrared monitor- 
ing may lead to reliable triggers for simultaneous multiwavelength campaigns to 
study jet formation processes. 

Subject headings: binaries (including multiple): close — stars: individual (4U 
1543-47) — black hole physics 



1. Introduction 

Soft X-ray transients (SXTs) are a subclass of X-ray binaries which comprise a compact 
object (black hole or neutron star) and a low-mass (< 1 Mq) secondary star. SXTs spend 
most of their time in quiescence but occasionally undergo X-ray outbursts during which 
the mass accretion rate through the disk increases dramatically and matter is accreted by 
the compact object. These objects provide an excellent means of studying accretion pro- 
cesses and how they relate to the fundamental binary parameters that can be measured in 
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quiescence. Jets have been observed in SXTs, particularly during low-hard states (Fender 
2001). SXTs are excellent candidates to answer outstanding questions about jet formation 
and behavior. To do so, however, one must obtain simultaneous (or quasi-simultaneous) 
multiwavelength data in order to relate the timing and origin of emission from various parts 
of the binary. 

The SXT, 4U 1543-47, was discovered on 17 August 1971 by the Uhuru satellite (Matil- 
sky et al. 1972). Since then additional X-ray outbursts have been observed in 1983, 1992 
and 2002. The optical counterpart was found during the 1983 outburst by Pederson et al. 
(1983). Dynamical studies of the companion star (Orosz ct al. 1998) show that the compact 
object is most likely a black hole with a mass of 2.7 < Mi < 7.5 M-q. 

The 2002 outburst began on MJD 52442 (UT June 16.683) as observed by RXTE 
(Miller and Remillard 2002; Park et al. 2004). During this outburst the soft X-ray light 
curve exhibited a classical fast rise, exponential decay profile as observed in many other 
SXTs (Chen, Shrader and Livio 1997). The top panel in Figure 1 shows the RXTE/ASM 
hght curve from MJD 52390 to 52550. At MJD 52443.5 the Molonglo Observatory Synthesis 
Telescope (MOST) detected a radio counterpart at 843 MHz with a flux density of 11 ± 1.5 
mJy (Hunstead and Webb 2002). From early on in the outburst until around MJD 52477, 
4U 1543-47 was in a high-soft state. A transition then occurred on this day into the low-hard 
state (Kalemci et al. 2002). Radio observations late in the outburst are also reported by 
Park et al. (2004). 

We observed the 2002 outburst of 4U 1543-47 as part of the YALO consortium (Bailyn 
et al. 1999)^. Our observations cover most of the outburst cycle in the infrared and some 
in the optical. A "secondary maximum" appeared in the optical and infrared (OIR) light 
curves which is not apparent in the RXTE/ASM data but, near its peak, corresponds to the 
second MOST radio detection. In this paper we present our OIR observations and reduction 
(§2), discuss the OIR light curve morphologies (§3) and focus our analysis on the secondary 
maximum (§4). In §5 we discuss our results and argue that the secondary maximum flux may 
be synchrotron emission which originates from a jet. We also note that secondary maxima 
may be a common phenomenon in SXTs during outburst decline after the source transitions 
into the low-hard state. In §6 we give our conclusions. 
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2. Observations, Data Reduction and Calibration 

V- and J-band images of 4U 1543-47 were obtained on a daily basis (weather permit- 
ting) from MJD 52323 (UT 2002 February 17) using the YALO 1.0m telescope together with 
ANDICAM^, a dual-channel imager capable of obtaining optical and infrared data simulta- 
neously. Data was recorded by a Lick/Loral-3 2048x2048 CCD on the optical channel and a 
Rockwell 1024x1024 HgCdTe "Hawaii" Array on the infrared channel. On UT 2002 June 5 
the YALO optical CCD failed. However, J-band observations continued on a daily basis. 

When the J-band hght curve was observed to have risen by ~ 0.5 mag (~ MJD 52440) 
we initiated additional observations in the K-hsaid. Optical observations were obtained on 
a daily basis, when possible, from UT 2002 June 18 to August 15 (MJD 52442.9 - 52500.8) 
using the 74 inch telescope at Mount Stromlo Observatory with the Cassegrain imager, a 2K 
X 4K CCD and B, V and I filters. Observations ceased when the object was no longer high 
enough in the night sky to observe. 

Exposure times for the optical frames were 300 sec for YALO images and 120 sec for 
Stromlo images. For the infrared, 7 images of 90 sec each were obtained in i^-band and 5 
images of 60 sec in J-band. Each image was shifted via an internal mirror in right ascension 
or declination by ~40 arcsec. 

Optical data reduction, including bias subtraction and flat fielding, was performed using 
the usual routines in IRAF. Photometry was performed using DAOPHOT in IRAF. Calibration 
was derived by calculating the offsets from a secondary star, the magnitude of which was 
measured previously (J. Orosz, private communication). The magnitudes were corrected for 
airmass extinction using coefficients taken from the CTIO table in IRAF. 

Infrared images were reduced using an in-house IRAF script which flat fields (using dome 
fiats), subtracts scaled sky images, shifts the images to a reference image then combines 
all images by averaging them. Aperture photometry was performed using qphot in IRAF. 
Calibration was done using the primary standard P9187^, observed on UT 2002 September 
2. Magnitudes were corrected for airmass extinction using atmospheric extinction coefficients 
taken from Frogel (1998). 
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3. Light Curve Morphology 

Figure 1 shows B—, V~, I—, J— and K-hand light curves together with the RXTE/ASM 
hght curve. Due to the lack of data in the optical bands it is more difficult to make qualitative 
or quantitative comparisons to the X-ray light curve. Therefore, most of our analysis of the 
outburst rise and peak deals with the IR bands. 

3.1. Start of Outburst 

The J-band light curve rises well before the RXTE/ASM flux becomes observable. In 
Figure 2 we have plotted the RXTE/ASM data (top panel) and J-band data (bottom panel) 
for the period leading up to, and just after, the start of the outburst. The mean value of 
the quiescent data was found to be 15.13 mag for J-band and counts/sec for the ASM. 
These are indicated by the dashed lines shown in Figure 2. To determine the start of the 
outburst we fitted the outburst part of the light curves with a hnear least-squares fit. For 
the X-rays we fitted data between MJD 52441.5 - 52444.7 (four points). For the J-band 
outburst data, two fits were made. The first included data from MJD 52435.5 - 52444.5, the 
second fit (shown as the solid line in Figure 2) was performed on data from MJD 52437.6 - 
52444.5. Including the two points at MJD 52435.6 and 52436.6 had a significant impact on 
the gradient of the line. This was the dominant error in determining the time of the start of 
the outburst. 

The J-band outburst began at MJD 52436.78 ± 0.58 (including the error mentioned 
above) whereas the X-ray outburst began at MJD 52441.51 ± 0.01. Fitting the X-ray data in 
log space resulted in an outburst start time of MJD 52439.38. Hence, the start of the X-ray 
outburst was delayed by ~ 3-5 days with respect to the infrared. Similar delays have been 
seen in other SXTs (Orosz et al. 1997). The slopes of the rise are 91.91 ± 0.01 counts/sec/day 
for the X-rays and -0.24 ± 0.01 mag/day for the J-band. Since the IR began to rise before 
the ASM, this suggests that the outburst was triggered in the outer regions of the accretion 
disk. Caution must be given, however, as the sensitivity level of the RXTE/ASM is orders 
of magnitude above the quiescent flux level for this source, so fine details of the X-ray light 
curve in the early stages of the outburst are undetectable. 

3.2. Outburst Maxima and Flares 

In comparing the J-band light curve to that of the ASM it is not apparent that there 
is any noticeable difference in the peak times. We determined the peak time in the J-band 
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light curve by fitting Gaussian profiles to the rising parts of the light curve. We fixed the 

?/-axis intercept in the fits to the mean values determined for the rise start times (§3.1). 
This was added as a constant to the Gaussian function. We restricted the fit between MJD 
52400-52450. The peak, amplitude and Gaussian width were varied. The best fit gave a 
peak time of MJD = 52447.1, amplitude = -1.9 mag and width = 4.8 days. The best fit 
Gaussian is plotted in Figure 3. A fit was not performed on the X-ray data as the small 
number of datapoints gave large errors. 

After the primary X-ray peak, the X-rays and OIR decline until ~ MJD 52460 when the 
K-hand light curve, and to a lesser extent the J-band, experiences a small fiare. Just before 
this fiare, the source transitioned briefiy into a steep power-law state and low-frequency 
QPOs were observed by Park et al. (2004). No optical or infrared data were obtained by 
us during the time of transition. At ~ MJD 52462, the K-hand returned to the previous 
decline rate from the primary outburst. Since the overall decline was not interrupted, we 
suspect that the small flare emission is due to a mechanism different from that of the primary 
outburst. 

At ~ MJD 52484 the J- and K-hand fluxes experience a signiflcant rise, which we call 
the "secondary maximum" . This rise in flux also appears in the optical bands but is much 
weaker. The peak strength of the secondary maximum increases with wavelength, that is, it 
is smallest in 5-band and highest in K-band. 

In Figure 4 we have plotted a zoomed part of the J- and K-band hght curves around the 
secondary maximum. One can separate the two components: a steadily declining flux from 
the primary outburst on top of which is superimposed additional flux from the secondary 
maximum. The first component we will call the "underlying fiux" . The secondary maximum 
is fiux additional and separate to that of the linearly-decaying fiux. Hence, it seems likely 
that its origin is also separate from the processes governing the primary outburst. 

We estimated the peak times in each band by fitting two-sided Gaussians to the sec- 
ondary maximum. Two-sided Gaussians were used as they provided a much better fit than a 
single Gaussian. The data were first detrended for the underlying fiux by removing a linear 
component fitted to the data either side of the secondary maximum. The best fit Gaussian 
parameters are shown in Table 1. The secondary maximum peaked in both wavebands at 
the same time: in the J-band, the fitted peak time was MJD 52486.31^1 (1-cr) while for 
K-hand it was MJD 52486.51?^ (l-a). 

We reiterate here that the second MOST radio detection was made near the peak of 
the secondary maximum. It is not known whether a radio source was present before or after 
this detection as observations were not conducted. Although there is no obvious change in 
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the RXTE/ASM light curve, there may be effects seen in the RXTE/PCA data which will 
be presented elsewhere (Kalemci et al., in prep.). 

After the secondary maximum, the OIR light curves return to the original decline rate 
of the primary maximum until reaching previous quiescent levels. 

4. Spectral Energy Distributions 

To understand the source of the secondary maximum emission, we have constructed a 
spectral energy distribution (SED) of the flux from the secondary maximum only at the time 
when the second MOST radio detection was made (MJD 52487). The OIR magnitudes of 4U 
1543-47 were first corrected for interstellar extinction using E{B — V)— 0.5 ± 0.05 (Orosz et 
al. 1998), Ay = 3.2 E{B - V) ^ 1.6 (Zombeck 1990), A^/Ay = 1.324 and Aj/Ay = 0.482, 
Aj/Ay = 0.282 and Ak/Av = 0.112 (Rickc and Lebofsky 1985). As there were no B, V or 
/ observations made on this date, wc extrapolated the optical light curves to estimate their 
magnitudes. For this date we estimate the secondary maximum magnitudes at the peak to 
be i?o=14.2, V^=14.2 and /o=13.9. Wc estimated the magnitude of the underlying flux at 
the time of the secondary maximum peak to be Bo = 14.5, Vo = 14.6 and lo = 14.6. In the 
IR bands, the secondary maximum peak is at Jo — 13.1 and Ko — 12.1 and the underlying 
flux magnitudes are Jo — 14.3 and Ko — 14.1. All magnitudes were converted to fluxes 
using zeropoint fluxes given by Allen (2000) and Campins, Rieke, and Lebofsky (1985). In 
each band, the underlying flux was subtracted from the peak flux giving the residual flux, 
corresponding to the flux of the secondary maximum only. 

The SED of the secondary maximum flux at its peak is presented in Figure 5. The top 
panel shows the data in log F,^ v. log v space and the bottom panel in log vF^, v. log v space. 
In the top panel the optical and IR points lie almost on a straight line. There is clearly a 
break somewhere between the K-hand point and the radio but it is impossible to say where 
exactly. In the bottom panel we can see that most of the power is in the OIR bands. 

4.1. Blackbody fits 

We attempted to flt the SED with a single blackbody, peaking at or around the K-haad 

point since this has the greatest flux. Two example fits arc shown in Figure 6. The solid line 
for = 1350°K peaks at K, while the dashed line is the best fit to the I, J and K data, with 
a blackbody temperature of T55 = 1650°K. Neither of these models fit the optical data. For 
a temperature of 1650°K, the emitting region would need to have a radius of ~ 1.4 x 10^"^ cm 
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given a distance of 9.1 kpc (Orosz et al. 1998). For an orbital period of 1.123 days (Orosz et 
al. 1998), the size of the primary Roche lobe is ~ 5 x lO^^cm. Hence this immediately rules 
out an optically thick accretion disk as the source of the observed flare radiation. 

The SED was fitted with multicolor-blackbody models (see Mitsuda et al. 1984, for 
parameter definitions). We fixed the inner radius, rj„ = Sr^, corresponding to the innermost 

stable orbit of a black hole (where ig is the gravitational radius). The mass of the compact 
object (black hole) was fixed at 5 Mq (Orosz ct al. 1998). In Figure 7 we show two fits: 
the first (solid line) is for a disk with Tout = lO'^'"^ 1^ (= 4.9 x lO^^cm) and M = 10^^ 
gs~^, the second (dashed line) is for Tout = 10^'^ rj„ (= 1.1 x lO^'^^cm) and M = 10^'^ gs~^. 
Although the models fit the /-, J- and i^T-band points, they cannot account for the B and 
V points. In the first case. Tout ~ 100 times the primary Roche lobe. In the second case. 
Tout is within the primary Roche lobe but M is far less than the average observed accretion 
rates of X-ray binaries, that is ~ 10^^-10^^ gs~^ (Frank, King, and Raine 2002). Hence 
multicolor-blackbody models cannot account for the secondary maximum flux. 

We fltted multicolor-blackbody models to the underlying outburst flux to check that it 
could be explained by an optically thick accretion disk. We chose to examine data at two 
times: MJD 52454 and MJD 52487. Data from the former time was taken from when the 
primary outburst had passed its peak but still bright relative to quiescent levels. It was also 
at this time that wc had data in all five bands. The latter time is that near the peak of the 
secondary maximum fiux. Wc constructed the SEDs following the method described here, 
subtracting the quiescent fiuxcs of 4U 1543-47 from the outburst fiuxcs thereby isolating 
the fiux related to the outburst source. Again, fiuxes were converted to mJy using the 
same zeropoint fluxes mentioned above. Errors for the SED points for MJD 52454 were 
calculated from photometric errors. Errors for the SED of the underlying flux at MJD 52487 
were estimated from the error in estimating the baseline connecting points either side of the 
secondary maximum. This was found to be ± 0.1 mag. 

When constructing multicolor-blackbody models we again flxed rj„ = Sr^ and Mi = 5 
M0. Figure 8 shows the SED at MJD 52454 and three typical multicolor-blackbody fits. 
We first chose ro„t = 10^rj„ which was within the primary Roche-lobe radius (~ 5 x 10^^ 
cm ~ 10^rj„). The resulting fit (shown as solid line) gave M = 10^^'^ gs~^ and Tout = 
5970°K. This M is somewhat higher than usually observed (see above), hence we chose M 
= 10^^'° gs~^ (dashed line) and 10^'''° gs^-*^ (dot-dashed line) for the next two fits. The 
corresponding temperatures of the outer disk were 5632°K and 5314°K, and Tout = 10^'^rj„ 
and iC'^rj^, respectively. In all three cases the models could provide a good overall match 
to the SED, although the values were never satisfactory, perhaps due to systematic errors 
from subtracting the flux due to the secondary star. All temperatures are very reasonable 
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given that the peak of the outburst has passed and that the accretion disk has probably 
transitioned to a cooler state at this time. The corresponding M are also plausible. 

The same fitting procedure was performed for the SED of the underlying flux at MJD 
52487, using the same three ro^t as above. The SED and resulting fits are shown in Figure 

9. Again we see good fits with the corresponding temperatures being 1125''K (solid line), 
1781''K (dashed hne) and IGSCK (dot-dashed line) and M = 10^^ gs"\ 10^^ gs"^ and 10^^ 
gs~^, respectively. Although these outer-disk temperatures seem low, fits to ellipsoidal light 
curves of other systems suggest that Tout can be < 2000°K (e.g. Beer and Podsiadlowski 
2002). 

4.2. Non-thermal emission 

Since the single-blackbody and multicolor-blackbody models cannot adequately describe 
the SED of the secondary maximum flux, we briefly look at non-thermal emission models. 

We fitted the SED of the secondary maximum fiux with a broken power-law in which 
each piece of the power-law is described by F,^ cx u^'" . We use this model to represent non- 
thermal emission from a synchrotron source. The SED and resulting fit are shown in Figure 

10. The solid line is the fit to the optical and IR points for which a = -0.79 ± 0.12. This 
power-law fits the optical/lR data well, in contrast to the thermal models. The dashed line 
is a power-law fit passing through the radio point assuming that the break is coincident with 
the i^-band point. In this case, a — 0.063 ± 0.003. This must be taken as a lower limit as 
it is possible that the turnover is at longer wavelengths. The broken power-law provides a 
much better fit than the thermal, blackbody models. Hence, we suggest that the source of 
the secondary maximum emission is synchrotron radiation. 

5. Discussion 

SXTs for which secondary maxima have been observed include A0620-00, GS 1124-68 
and GRO J0422+32 (Chen, Shrader and Livio 1997) and XTE J1550-564 (Jain et al. 2001). 
In the first three objects, the maxima were observed at short wavelengths {B and UV) and 
were quite weak. The secondary maximum appeared ~ 40-80 days after the primary X-ray 
maximum. The simultaneous soft X-ray light curves, which had not yet reached quiescent 
levels at the time, also showed secondary maxima. The soft X-ray light curve of 4U 1543-47 
in 1971 (Li, Sprott and Clark 1976) exhibited a secondary maximum ~ 70 days after the 
primary peak. No optical nor infrared observations were made at the time. 
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XTE J1550-564 (Jain et al. 2001) exhibited very similar behavior to that of 4U 1543-47 
presented in this paper. The secondary maximum appeared ~ 50 days after the the X-ray 
maximum and increased in strength from optical ( V- and /-band) to infrared (//-band) 
wavelengths. The ASM had already reached quiescent levels and did not show a flare at the 
time of the secondary maximum. There was, however, a deviation in the exponential decay 
of the RXTE pointed observations during the secondary maximum (Tomsick, Corbel and 
Kaaret 2001). Corbel et al. (2001) reported a radio detection near the peak of the secondary 
maximum. The detection was made in four bands giving a spectral index of 0.37 ± 0.10. 

The secondary maxima in these sources have characteristics consistent with the one we 
have observed in 4U 1543-47. In particular, the strength of the peak increases from the 
optical to infrared wavelengths. The times at which the secondary maximum appears after 
the primary X-ray peak agree within an order of magnitude. This timescale could relate to 
some fundamental parameter of the binary system, such as the orbital period or mass of the 
compact object. It seems likely, however, that the trigger of the secondary maxima is related 
to the transition into the low-hard state. The secondary maxima of 4U 1543-47 and XTE 
J 1550-564 begin to rise ~ 7 and 8 days, respectively, after the transition into the low-hard 
state. Both objects exhibited radio emission associated with a jet which was detected near 
the peak of the secondary maxima. Jets are believed to form during the low-hard state 
(Fender 2001). If the secondary maximum emission does in fact originate from jets then one 
would expect to see it after the object has transitioned into the low-hard state. 



6. Conclusion 

We have presented OIR light curves of the 2002 outburst of 4U 1543-47. A strong 
IR flare appears soon after the source transitions into the low-hard state. Blackbody and 
broken power-law fits suggest that synchrotron emission is the most likely source of the 
secondary maximum fiux. The detection of a radio source at the peak of this flare supports 
this conclusion. Indeed, the presence of jets during this X-ray state would not be surprising 
given that we often observe jets in this accretion phase (Fender 2001). 

Since the secondary maximum is so prominent in the IR, it is clear that IR monitor- 
ing of SXTs can successfully trigger other simultaneous observations at other wavelengths 
and gather useful information into the formation of and flux contribution of jets. Future 
observations of such phenomena can be further improved by obtaining simultaneous radio 
observations in two bands so that a confldent flt of the power-law gradient may be obtained. 
UV observations will enable us to ascertain any upturn in the SED due to the dominance of 
an accretion disk blackbody component, although observations in this band are quite difficult 
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due to strong interstellar extinction. X-ray spectra would also provide valuable information 
on the relationship to the OIR data and if there is any contribution from thermal IC of disk 
photons. 
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Fig. 1. — Optical and infrared data of 4U 1543-47 in 2002. Tlic top panel is the one-day averaged 
RXTE ASM data provided by the ASM/RXTE team. 
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Fig. 2. — Quiescent period and early part of the 2002 X-ray outburst of 4U 1543-47. The 
top panel is the one-day averaged RXTE/ASM data, bottom panel is J-band data. The 
dashed line is the mean quiescent value. The sohd line is the fit to the outburst data. The 
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Fig. 3. — Gaussian fit to the initial rise in J-band light curve. The vertical-dashed line 
indicates best fit peak time of MJD 52477. The solid fine is the rising portion of best-fit 
Gaussian. 
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Fig. 4. — Two-sided Gaussian fits to the secondary maximium in J- and K-bands. Dashed 
hne is best-fit slope of underlying flux, sofid line is best-fit two-sided Gaussian. 
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Fig. 5. — Spectral energy distribution of the secondary maximum flux on MJD 52487 (UT 
2002 July 31). The B, V and / points were extrapolated from their respective hght curves. 
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Fig. 6. — Single-temperature blackbody fits to the optical and infrared section of the SED of 
the secondary maximum flux. Dashed hne is for T^b = 1650°K, solid line for T^^ = 1350°K. 
These models cannot account for the optical/infrared flux observed. 
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Fig. 7. — Two examples of multicolor-blackbody model fits to the SED of the secondary 
maximum flux. See text for model parameters. The multicolor-blackbody model is slightly 
better than the single-temperature blackbody model (Fig. 6) but there is still a large excess 
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Fig. 8. — Multicolor-blackbody model fits to the SED at MJD 52454 near the primary 
maximum. See text for model parameters. These models fit the data more satisfactorily 
than that seen for the SED of the secondary maximum flux (Fig. 7) and give reasonable 




Fig. 9. — Multicolor-blackbody model fits to the SED of tlie underlying flux of the secondary 
maximum at MJD 52487. As for the SED at MJD 52454, these models describe the data 
well. See text for model parameters. 
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Fig. 10. — Broken power-law fit to the SED. The sohd hne passing through the optical and infrared 
points has a slope of -0.79. The slope clearly changes between the radio and if -band points but 
we can not say where exactly. The dashed line joining the radio and if-band points has a slope of 
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Table 1: Two-sided Gaussian parameters for best fits to secondary maxima in J- and iiT-band. 
Errors shown are 1-a. 



Waveband 


Amplitude 

(mag) 


Peak time 
(MJD) 


(days) 


0-2 

(days) 


J 


-i.lU_0.20 


52486.3tiJ 






K 


-2 1 n+°-22 

'^■-'-'-'-0.23 


52486.51^;^ 


9 9+1.6 
^■^-1.6 


6.611^ 



